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TWO-POINT CORRELATION FUNCTIONS AND UNIVERSALITY FOR THE ZEROS OF SYSTEMS 
OF SO(N+l)-INVARIANT GAUSSIAN RANDOM POLYNOMIALS 

PAVEL M. BLEHERi, YUSHIAND ROLAND K. W. ROEDERi 


Abstract. We study the two-point congelation functions for the zeroes of systems of SO{n + l)-invariant Gaussian 
random polynomials on and systems of Isom(R^)-invariant Gaussian analytic functions. Our result reflects the same 
“repelling,” “neutral,” and “attracting” short-distance asymptotic behavior, depending on the dimension, as was discovered 
in the complex case by Bleher, Shiftman, and Zelditch. 

We then prove that the correlation function for the Isom(R^)-invariant Gaussian analytic functions is “universal,” de¬ 
scribing the scaling limit of the correlation function for the restriction of systems of the SO{k + l)-invariant Gaussian 
random polynomials to any n-dimensional submanifold M C RF^. This provides a real counterpart to the universality 
results that were proved in the complex case by Bleher, Shiffman, and Zelditch. 


1. Introduction 


This paper concerns the 50(ri+l)-invariant ensemble of Gaussian random polynomials on and the Isom(E^)- 

invariant ensemble of Gaussian random analytic functions on The SO{n + l)-invariant ensemble consists of 
random polynomials of the form: 

( 1 ) E 

where X € R"+^ and the Ua are independent and identically distributed (iid) on the standard normal distribution, 
Af (0,1). Here, we use the following multi-index notation: for any a G (Z>o)"^^, one defines: 


( 2 ) 




n+1 


2=1 


n+1 

|q:| = ai and 

i=l 



d\ 


n+1 


n aj\ 


We will study the simultaneous zeroes on the projective space KP" of the systems: 

(3) where F = (Fi(X), F 2 (X),..., F„(X)), 


where each Fi is an independently chosen random function of the form in Equation (1). Almost surely, the common 
zero set of F will be finitely many points. We equip RP" with the Riemannian metric obtained from its double 
cover by the unit sphere S" C R"+^. The simultaneous zeroes of ensemble (3) are invariant under the isometries 
by elements of SO(n + 1); see Section 2. Because of this symmetry, authors have described this ensemble as the 
“most natural” ensemble of a random polynomials defined on RP". For this reason, it has been extensively studied by 
Kostlan-Edelman[12], Shub-Smale[32], and others. 

The Isom(R")-invariant ensemble of Gaussian random analytic functions is defined by the following: 


(4) 


/ : R" -)> R" where / = (/i(a;), f 2 {x),fn{x )), 


with f,{x) := y] 




where Oq, are iid on the standard normal distribution, A/^(0,1). We will show in Section 2 that the zeroes of this 
ensemble are invariant under all isometries of R". We will see shortly that this ensemble is intimately tied to the 
SOin + l)-invariance ensemble in the scaling limit as the degree d —>■ oo. 
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The probability density of the zeros of the system (3) aix G is defined to be 
(5) p(x) = lim (3 a zero of F in iV,(x)), 

where Ns{x) := {y G K.P" : dist(a:, y) < 5}. It follows from the invariance that this ensemble (3) has a constant 
density of zeroes given by 


( 6 ) 


Pd{x) = TT 


n + 1 




see, for example, [12, Sec. 7.2]. Note: the volume of the real projective space is Tr^~T so the expected 

number of zeroes is simply . The analogous definition applies to the ensemble (4) which, because of the invariance 
under isometries of R" has constant density 


(7) 


p{x) = 




The correlation function between the zeros of the system (3) at the two points x and y in RP" is defined to be 


( 8 ) 


Kn,d{x,y) := lim 


Pr (3 a zero of F in Ns{x) and 3 a zero of F in Ns{y)) 


s^o Pr (3 a zero of F in Ns{x)) Pr (3 a zero of F in Ns{y)) 


It follows from the SO{n + 1) invariance that Kn,d{x, y) depends only on the distance between x and y. For this 
reason, we can write Kn,d{x,y) = Kn,d{t), where t = distBp" (a;, y). Similarly, for any x,y G R", the two point 
correlation function ICn{x, y) between zeros of (4) depends only on distR" {x,y). We have 


Theorem 1. For any x ^ y G RP", let t = distBpzi [x, y). For fixed d > 3, the correlation function between zeros of 
the SO(n + l)-invariant ensemble satisfies 


(9) Kn4{x,y) = Kn,d{t) = + ”) as t-t 0, where A„,d 


d-i \ V^r(^) 

dt J 2r(^) 


Theorem 2. For any x y G R", let t = distRn (a;, y). The correlation function between zeros of the Isom(R"')- 
invariant ensemble satisfies 


( 10 ) 


Kn{x,y) = Kn{t) = An ” + O ") as f —>• 0, where An 


2r(^) 


and 

(11) Kn{t) = 1 + O ast^oo. 


Given a submanifold M C RP^ having dimension n, the restrictions of n of the polynomials chosen iid from 
the SOik + l)-invariant ensemble (1) has a well-defined zero set which again consists a.s. of finitely many points. We 
give M C RP^ the metric induced by the double cover of RP^ by the unit sphere S^. More specifically, we obtain 
a Riemannian metric on M using the inclusion of tangent spaces TpM C TpRP^. When restricted to a sufficiently 
small neighborhood of the origin, the orthogonal projection projp : TpM -G- M provides a system of local coordinates 
on M. We will use these systems of local coordinates to study the correlation between zeros of the restriction of the 
SO{k -b 1) invariant ensemble to M. The next theorem expresses Kn{x, y) as the universal correlation function in 
the scaling limit d —>■ oo for the restriction of the SO{k -b l)-invariant ensemble to M C RP^: 


Theorem 3. Let M C RP^ be a (7^ submanifold of dimension n and Kn,d,M{x, y) denote the correlation function 
between zeros of n polynomials chosen iid from the degree d SO{k + 1) invariant ensemble restricted to M. Then, for 
any p G M and any x,y G TpM we have 


Kn 


,d,M I projp 


=Kn{^^y) + o(4i 


Vd 


\VdJ 


s/d 


d —^ 


The constant in the estimate is uniform on compact subsets of TpM x TpM \ Diag, where Diag = {(x,y) G 
TpM X TpM : x = y}. 
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Our techniques are largely based on those of Bleher and Di [4], who use the Kac-Rice formula (see Section 3 
below) to study the n-point correlation functions for the iS'0(l, 1) and 5'0(2)-invariant polynomials in one variable. 
Moreover, our results in the higher dimensional real case yield the exact same short-distance asymptotic behavior 
(with a different constant) as those of Bleher, Shiffman, and Zelditch [7, 8, 9] in the complex case. These asymptotic 
behaviors can be interpreted as “repelling” for n = 1, “neutral” for n = 2, and “attracting” for n > 3. See Figure 1 
for numerical plots of Kn{t) for n = 1, n = 2, and n = 3. 

We remark that calculation of the leading order asymptotics is more delicate in the real case than in the complex 
case because one cannot apply Wick’s Theorem to the real Kac-Rice formula. Similar types of analysis have been 
done in the real setting by Nicolaescu, who studied the critical points for random Fourier series. It is interesting that 
he found the same exponent of 2 — n arising in his work [29, Eqn. 1.15]. 

Theorem 3 above provides a real analog of the universality results that were obtained in the complex setting by 
Bleher, Shiffman, and Zelditch [7, 8]. Thus, the plots shown in Figure 1 depict the universal scaling limits of the 
correlation functions for any submanifold M C RP^ of dimension 1, 2, or 3. 

The scaling limit used in Theorem 3 is needed to get a universal correlation function. This is illustrated in Section 8 
where we show that when restricted to a parabola y = hx^ the leading term from the correlation between zeros for the 
S'0(3)-invariant polynomials of degree 3 near a; = 0 depends non-trivially on b. More generally, it can be interesting 
to ask how the geometry of M affects the correlation function K for finite degree d. 


The proof of Theorem 3 easily adapts to to complex setting: The SU{k + l)-invariant ensemble of polynomials are 
obtained by interpreting the variables in (1) as complex and replacing the real Gaussians Ua with complex Gaussians. 
The Isom(C”)-invariant ensemble of Gaussian analytic functions on C"* is obtained by making the same adaptations 
to (4). We obtain: 


Theorem 4. Let M C CP^ be an complex analytic submanifold of dimension n and Kn^d,M{x^ y) denote the corre¬ 
lation function between zeros of n polynomials chosen iidfrom the degree d SU (fc -f 1) invariant ensemble restricted 
to M. Then, for any p £ M and any x,y £ TpM we have 


d^n,d,M 




Kn{x,y) -f o 



as d ^ oo. 


The constant in the estimate is uniform on compact subsets of TpM x TpM \ Diag, where Diag = {{x,y) £ 
TpM X TpM : x = y}. 


This serves as a weaker version of the results from [7, 8] in that M is required to be embedded in projective space 
(instead of being an arbitrary Kahler manifold), the line bundle is the hyperplane bundle (corresponding to the 
SU{n -f 1)-invariant ensemble), and only two-point correlation functions are considered. On the other hand, in the 
work of [7, 8] the manifold M is assumed to be compact. No such assumption is made in Theorems 3 and 4. For 
example, they can be applied at any smooth point p of a singular projective variety. 


For general background on Gaussian random analytic functions and polynomials, we refer the reader to [2, 3, 19, 
20, 32] and their references therein. Specifically to correlation functions, we refer the reader to the three papers listed 
above in the previous paragraph, as well as the works of Bogomolny, Bohigas, and Leboeuf [10], Tao and Vu [33], 
Bleher and Ridzal [6], and Bleher and Di [5]. 

Our work fits in within the context of the emerging field “random real algebraic geometry.” For example. Theorem 3 
applies to the restriction of the SO{k -f 1) ensemble to the smooth locus of a real-algebraic subset of RP^. We 
refer the reader to the works of Kostlan [24], Shub-Smale [32], Ibragimov-Podkorytov [21], Burgisser [11], Gayet- 
Welschinger [14], Ibragimov-Zaporozhets [22], Nastasescu [28], Gayet-Welschinger [15, 16], Lerario-Lundberg [25], 
Gayet-Welschinger [17, 18], and Fyodorov-Lerario-Lundberg [13]. 

The remainder of the paper will be organized as follows: In the following Section 2, we study the invariance 
properties of the ensembles from (3) and (4). We then use the invariance to reduce Theorems 1 and 2 to suitable 
versions in affine coordinates (Theorem 9). In Section 3, we recall the Kac-Rice Formulae for the density and for the 
correlation functions, the main tools used in our proof. In Section 4 we compute the covariance matrices needed to 
prove Theorem 9, as well as their determinants, inverses, etc. Theorem 9 consists of two statements (short-distance 
asymptotics and long-distance asymptotics), which are proved in Sections 5 and 6 respectively. Section 7 is dedicated 
to proving Theorem 3 about universality of the scaling limit. Section 8 provides an example showing that for finite 
degree the leading asymptotics depends on the geometry of the submanifold M C RP^. In Section 9 we explain the 
changes that need to be made to the proof of Theorem 3 in order to prove the complex version. Theorem 4. 
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n=1 



n=2 



n=3 



Figure 1. Universal two-point limiting correlation functions Kn {t) for n = 1,2, and 3, demon¬ 
strating the repelling, neutral and attracting behaviors. For n = 1, the graph is obtained from 
Formula (5.35) in [4]. For n = 2 and n = 3, the graphs were computed using Monte Carlo in¬ 
tegration applied to formula (86) with 10^ and 10® points, respectively, for each t. The data was 
smoothed out by replacing each value with the average of it and the 14 nearest neighboring points. 
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Appendix A contains the proof of a general estimate which is used in Sections 6 and 7. In Appendix B, we prove a 
result regarding the volume of random parallelotopes which is needed in Section 5. 

Notation: Let diag^. (A) denote the block-diagonal matrix with k copies of the square matrix A along the diagonal. 


2. Invariance Properties and Reduction of Theorems 1 and 2 to local coordinates 

The SO{n + l)-invariant ensemble and the Isom(]R”)-invariant ensemble are instances of the following definition; 

Definition 5. A Gaussian analytic function h : R" — >■ R™ is an m-tuple {hi{x), h 2 {x ),..., hm{x)) of functions 
hi : R" —>■ R chosen iid of the form 

(12) hi{x^ .— ^ ^ CgOaX , 

Ct 

where the Ua are chosen iid on the standard normal distribution A/" (0,1) and the coefficients Ca are chosen so that 
^ ColX°^ converges for all x £ R". 

a. 

Lemma 6. A Gaussian analytic function h : R" — >■ R"* almost surely converges uniformly on compact subsets o/R" 
and moreover is real analytic on R". 

Proof The proof of Lemma 2.2.3 from [20] applies to show (12) almost surely converges uniformly on compact 
subsets of C" and hence defines a random complex analytic function on C". By restricting the resulting functions to 
R", we obtain the desired result. □ 

In particular, Lemma 6 justifies our consideration the Isom(R")-invariant ensemble (4) as actually defining a ran¬ 
dom function. The following two lemmas justify our terminology “SO{n + 1) invariant ensemble” and ‘Tsom(R”)- 
invariant ensemble”: 


Lemma 7. The zeroes of the system F given in (3) are invariant under the action of SO{n + 1). That is, for any open 
set U C RP" and any A £ SO{n + 1), we have Pr {F has a zero in U) = Pr {F has a zero in A (U)). 

Proof Each Fi{X) defines a Gaussian process on R”+^, with mean 0 and covariance function 
(13) E{F,{X)F,{Y)) = ■ ^)''- 

Since any Gaussian process is uniquely determined by its first and second moments [19, Theorem 2.1], this process is 
invariant under SO{n -f 1). Therefore, the zeros within RP" are also invariant under the action of SO{n -I- 1). □ 

Propositions. The zeroes of the system f = {fi,f^) from (4) are invariant under any isometry of MA. That is, for 
any open set U C R" and any isometry I : R" —> R", we have Pr (/ has a zero in U) = Pr (/ has a zero in I (U)). 


Proof The zeroes of f are the same as those of 

(14) g := {gi,g 2 ,... ,gn) where g,{x) := f,{x). 

Each gi{x) defines a Gaussian process on R", with mean 0 and covariance function 

n / OO / 


(15) E {g,{x)g,{y)) = e ^ 


- ii=>^ir+iiyii 


E 




— Q-2\\^-y\\ 


a! 




The result follows because (15) is clearly invariant under isometries of R". 


□ 


We will now use these invariance properties to reduce the proofs of Theorems 1 and 2 to a particularly simple pairs 
of points and to local coordinates. The two points 


(16) 


X = 


1 : 0 : 


t 


; 0 : -- 


and y = 


1 : 0 ; ••• : 0 : 


t 


(given here in homogeneous coordinates) have distance 

distRp" {x,y) = 2 arctan 


— t Oif^') as t —y 0. 


Thus, in order to prove Theorem 1, it suffices to verify (9) for this pair of points. 
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Note that (xi, ..., Xn) ^ [1 : xi : : Xn\ provides a system of local coordinates in a neighborhood of x and y. 

In these coordinates, the SO{n + l)-invariant ensemble becomes 

fd = {fdAx)Jd, 2 ix), fd,nA)) 1 

where each ^ is chosen independently of the form 

(17) /.(*) = ^ 

|a|<d V \ / 

i=l 

and the are iid on the standard normal distribution V (0,1). 

In summary: Let /C„ d(a:, 2 /) and ICn{x,y) denote the correlation functions between zeros of the SO{n + 1)- 
invariant ensemble, expressed in affine coordinates (17), and between zeros of the Isom(K")-invariant ensemble (4), 
respectively, and let 

(18) :=/C„.d((0,...,0,-V2),(0,...,0,f/2)) and /C„(f) :=/C„((0,..., 0,-f/2), (0,..., 0, f/2)) 

In order to prove Theorems 1 and 2, it suffices to prove: 

Theorem 9. 

(1) As t ^ 0 we have 

ICn,d{t)=An,dt^-" + OA-A, and ICAt) = 

where An^d and An are given in (9) and (10), respectively. 

(2) As t ^ oo we have 

~ 1+0 ^ . 


dn X 


where 


d\ 


(d- |q:|)! n 


3. Kac-Rice Formula 


The main technique used in this paper is a variant of the classical Kac-Rice Formula [23, 30, 31 ] that was developed 
for correlations between zeros of multivariable Gaussian analytic functions by Bleher, Shiffman, and Zelditch in [8, 
Section 2]. 

We will begin this section with a statement and proof of the Kac-Rice formula for the m point correlation function 
with m arbitrary. At the end of the section we will rephrase the results as needed in this paper. Let h : R" —R" be 
a Gaussian analytic function and let ..., a:™ be m distinct points in R”. The m-point correlation function for the 
zeros of h is 

Pr (3 a zero of h in Ns(x'^) for each i = 1,..., m) 


(19) 


Kn(x ^,..., a;™) := lim 
5—>-0 


n™ 1 Pr (3 a zero of h. in iV 5 (a;®)) ’ 

where Ns{x'^) is the ball of radius <5 > 0 centered at a;® € R”. 

The m-point density for the zeros of h is 

(20) Knix^,..., a;"®) := lim -—-—Pr (3 a zero of h in Ns{x'^) for each z = 1,..., m) 

ny , , j (VsA ^ / 


/ V is the volume of each ball. When m = 1, Kn{x) = Pn{x), the probability density (5). For 
r(^-i-i) 

m > 1 we have 


where Vs = 


( 21 ) 


Kn{x^ 


,xn = 


1 


-Kn(x\ . . . ,XA- 


Consider the Gaussian random (mn^ 
(22) v:= Vh 2 (x^) ... 


nZl Pn(x^) 

TOn)-dimensional column vector 
Vhn(x^) . Vhi(x"^) Vh2(x"^) 


h(x^ 


h(x^) 


Vhn(x’^) 

h(x^)]\ 


where each gradient vector V/ii(a:^) and each vector h{x^) is concatenated into v in the indicated location. 
Let ..., be n X n matrices whose rows are ... ^^ for each z = 1,..., m. Let 


« = f 




tl rr. 

Sn Sn 


(23) 
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be the mn^ dimensional column vector formed by concatenating the first rows of each of the matrices ..., 
followed by their second rows, etc. 


2 , 

Proposition 10 (General Kac-Rice Formula). Suppose the covariance matrix C = of the vec¬ 

tor (22) is positive definite. Then, the m-point density Knix^, • ■ •, x™") for the zeroes of the system h is given by 


(24) 


Kix\...,x^) 


1 

(2^)'""(”+i)/" v/ditC 




where is the mv? x mvf principal minor ofC ^ and u is as in (23). 


Proposition 10 is easily obtained from [8, Theorem 2.2] by using the suitable Gaussian density Dk{0, z) in their 
formula (38). In order for this paper to be relatively self contained, we present a proof of Proposition 10 below. 

We start with the following lemma. Define the derivative Dh{x) as a linear map w i— Dw, where D is the matrix 


and w = [wi,..., WnY. 


D = 


dhj 

dxi 


(x) 


i,j=t 


Lemma 11. We have 

(25) Kn{x^ ,..., a;"*) = lim ■^^^-y^Pr^h,(a;®) € Dh{xY{Ns{0)) for each i = 1,..., . 


Proof. We begin by cutting off the tails of h. Let i? > 0 be chosen sufficiently large so that for all (5 > 0 sufficiently 
small and alH = 1,..., m we have Ns{x'^) C {||at|| < i?}. Consider the following bounds on the derivatives of h: 


(26) 


||Dft.(a:*)|| < A for all 1 < i < m and 


d"^hj{x) 

dxkdxi 


< A for all 1 < j,k,l < n and all ||a;|| < R. 


For any Gaussian analytic function h 

Pr (ji satisfies condition (26) ^ 1 as A — c». 

Therefore, it will be sufficient to prove (25) under the hypotheses (26). (The constant A > 0 will be fixed for the 
remainder of the proof.) 

2 

As in the statement of Proposition 10, let be the n x n matrices, for i = 1,... ,m, and let u G R™" be given 
as in (23). Let s := [ • • • s™ ] ^ be the mn-dimensional column vector, where each s® G R". For any 

open subset U C R" having compact closure and any e > 0 let 


:= {x G (7 : dist(a;, d(7) > e} and 17^ ■= G R" : dist(at, 17) < e}. 


For any i) > 0 and B > 0 consider 

(27) := G : ||C®|| < A and s® G (l7iV5(0)))±^, for each i = 1,..., m} and 

(28) := |(m,s) G : ||r|| < A and s® G ^®(A^5(0)) for each i = 1,..., m| . 

Here, ^'^{Ns{0)) denotes the image of the ball Ns{0) under the the linear map from R" to R" expressed in terms of 
the standard basis on R" by the n x n matrix 

The bounds ||^®|| < A, for i = 1,..., m, imply that as (5 —0 for fixed B > Owe have 

(29) Yo\{E+^\Ef^) = 0{S^-+n- 

Let 


Hs := {h : h has a zero in Ns{x’') for each i = 1,..., m} , 

■ {Dh{xY,...,Dh(^x'^),h{x^),...,h{x"^)) G Efg^ and 
H^:={h : {Dh{xY,.. ■, Dh{x^),h{xY, ■.. ,h{x^)) G E^} . 

It is immediate from the definition of the sets Efg and Eg that for any (5 > 0 and any B > 0 that 
(30) C C 
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Meanwhile, by assumption (26), Taylor’s Theorem gives that there exists B > 0 (depending on the bound A) so that 
for all sufficiently small i5 > 0 we have 

(31) cHsC H+^. 

Since (^Dh{x^),..., Dhix^), h{x^),..., is a Gaussian random vector, its probability distribution is ab¬ 

solutely continuous. Therefore, (29) implies that 

(32) Pr (/r € H+^ \ 

Combined with (30) and (31) this implies that under the assumption (26) we have 


(33) 


since 


lim -—-— Pr = lim -—-— Pr fft, £ Hf \, 

5^0 (Pa)™ V V V V’ 


(Vs)'^ is bounded below by a constant times b'"". 


□ 


Proof of Proposition 10. We use Lemma 11 to replace the definition of Kn with (25). Using the formula for 
Gaussian density, we have 

(34) 

^ ^ -41 c- 


Kn{x^, . . . ,X^) = 


lim 


1 


(27r)’"”("+i)/2 VditC ^^0 (Vs) 


u 

s 


u 

s 


dsdu. 


n2 «l(Ar5(0))x-XC™(AfP0)) 


where 


u 

s 


denotes the column vector obtained by stacking the two column vectors u and s. 

Because C is positive definite, the integrand decays rapidly at infinity. Thus, the Dominated Convergence Theorem 
allows us to interchange the limit with the first integral. We also multiply and divide by Y\aLi I ^*|, obtaining: 

(35) 


= 


lim 


nr=iidetei 


-4 c 


(2^)'"”(”+i)/ 2 Tditc J s^oUtiyoieiNsiO)) 


u 

s 


u 

s 


dsdu. 


^mn 


^i(iV5(0))x-x^-(iV5(0)) 


The Integral Mean Value Theorem implies that 
(36) 


-4 c 


lim - / 

^^oi\T=iyoie{Nsm J 

e{N5{0))x-xi^{Ns{0)) 

which completes the proof. 


u 

s 


u 

s 


-4 c 


ds = e 


u 

0 


u 

0 




□ 


Remark 12. Proposition 10 shows that the correlation measure is absolutely continuous off of the “diagonal" where 
x'’ = xl for some i j (hence the name “correlation function"). Thus, in the definition (19) of Kn{x ^,..., at"*) one 
need not use round balls Ns(x'^). Rather, any sequence of shrinking neighborhoods of each x* suitable for computing 
a Radon-Nikodym derivative will suffice. 

On certain occasions we will need the following lemma, which is proved in Appendix A, to make estimates involv¬ 
ing the Kac-Rice formula (24). 


Lemma 13. For any positive definite mnf x rrcnf matrix A 


(37) 


nidetrie-^ 






2=1 


= of||A-B||i/2 
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for any mn? X mn^ matrix B sufficiently close to A. (Here u is as in (23) and || ||oo denotes the maximum entry of 
the matrix.) 


Remark 14. After reading a preprint of this paper, L. Nicolaescu informed us that a very similar estimate appears in 
his paper [Prop. A.l ][29], whose proof was provided by G. Lowther in a discussion on Math Overflow [26]. 


We close the section with simplified rephrasings of Proposition 10 in the cases m = 1 and m = 2 that will be used 
throughout this paper. In both cases it will be better to order the random vector v from (22) in a way that will make the 
covariance matrix C block diagonal. (It will come at the cost of the definition for being slightly more complicated.) 
If TO = 1 we reorder the random vector v as 


(38) u := [ hi{x) V/ii(£c) ... hn{x) Vhn{x) ]^. 

2 I 

Because the components of h are chosen iid, the resulting covariance matrix C = {EviVjfl Aff will be of the form 
C = diag„(C), where C = {EviVj)'^A^i and diag„ denotes the block diagonal matrix with n copies of C on 

the diagonal. The vector u becomes it = [ ... ^ where ^ is a n x n matrix. 

Proposition 15. (Kac-Rice for Density) Suppose the covariance matrix C = {EviVj)'^^^^'^ of the vector (38) is 
positive definite. Then, the density of zeroes of the system h is: 


(39) 


VdeFC y 


where Jl is the matrix of the elements of C ^ left after removing the rows and columns with indices congruent to 1 
modulo n + 1. 


If TO = 2 we reorder v as 

(40) V := [ hi{x) Vfii(a:) hi{y) Vhi{y) ... hn{x) Vhn{x) hn{y) 

which will again make its covariance matrix block diagonal C = diag„(C). Let ^ and rj be the nx n matrices whose 
rows are ... ^„ and yi,... T 7 „, respectively. Let M = [ r/i ^2 '<12 ■■■ Vn ]'be the vector formed 

by alternating the vectors and ry. 

Proposition 16 (Two Point Kac-Rice Formula). Suppose the covariance matrix C = of the vec¬ 

tor (40) is positive definite. Then, the two-point correlation function for the zeroes of the system h is: 

where ft is the matrix of the elements of C~^ left after removing the rows and columns with indices congruent to 1 
modulo n + 1. 


4. Calculation of the covariance matrices, their inverses, and ft 

Let Cn^d = Cnffit) and C„ = C„(t) be the covariance matrix for vector (40) applied to fd (Equation 17) and f 
(Equation 4), respectively, at the points 


( 42 ) 31 = (^ 0 ,..., 0,-0 

Lemma 17. Both Cn,d ond C„ are of the form 


and y = 0,..., 0, ^ 

A+ BT 
B A ’ 


(43) 


C:= diagJC), 


with 


C = 
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where A± and B are the following (n + 1) x (n + 1) matrices: 




a 

0 .. 

. 0 

±5' 


’ F 

0 . 

. 0 

ly 



0 

/3 ■ 


0 


0 

11 


0 

(44) 

A± = 





, B = 







0 


■ /3 

0 


0 


• V 

0 




0 .. 

. 0 

7 . 


— V 

0 . 

. 0 

T 


and where a,/3, <5, 7 ,/i, 77 , i/, and r are functions of d and t expressed in (50-57) for Cn,d cmd the functions oft 
expressed in (62-69) for C„. 

Proof Since the coefficients of fd^ and fdj (respectively fi and fj) are independent when i ^ j, only the entries of 
Cn,d (respectively C„) with i = j will have nonzero values. Thus, the covariance matrices will have the following 
block-diagonal structure; 


(45) 


= diag„(Cn.d) and C„ = diag„(C„), 


where Cn.d corresponds to the first 2n + 2 entries of v (and similarly for Cn). These entries correspond to 1 and 
/i, respectively. For ease of notation, we’ll drop the subscript 1: fd = fd.i, f = fi- 
For any z,w G M” we have: 


(46) 


E{fd{z)fd{w)) = {l + z-wf 


(47) E fd{z) 


dfd{'w)\ dE{fd{z)fd{w)) 


dw,. 


dwi 


= d Zi{l-\- z ■ w) 


d-l 


(dfd{.z) dfd{w)\ d‘^E{fd{z)fd{w)) d-l , , ^d -2 ^ 

(48) E{ —- = d [1z ■ w) d[d — 1)ZiWi [1z ■ w) , and 


\ dzi dwi 


dzjdwi 


(dfd{z) dfd{w)\ d‘^E{fd{z)fd{w)) d -2 ^ . 

(49) E I —-- ) =--= d(d — l)zjWi [1 z ■ w) for I J. 


V dzi dwi 


dzidwj 


Recalling that x = (O,..., 0, — |) and y = (O,..., 0, |) we now use expressions (46-49) to compute that the only 
non-zero covariances in Cn,d are 


(50) a := E{fd{x)fd{x)) = E{fd{y)fd{y)) = ( 1 + ^ 


(51) 5-.= e[ fd{x) 


, dfd{x) 

dxn 




(52) H E y ) = E y ) =‘‘{' + 1 

(53) 

\ dXn OXn J \ ayn ayn J \ 4 / 


d-l 


i -1 


for i n, 


d-l 


-I- -d(d- l)f^ ( 1 -f — 


d-2 


= -d- 1 - 


(54) p:= E{fd{x)fdiy)) = {l-\-x-y) = - — 

(55) . =-£(/, 

This proves that Cn d has the structure stated in Lemma 17. 


t / .2X d-l 


for i n, and 


1 - 


d-2 
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For Isom(K")-invariant ensemble f and any z,w G K" we have: 


(58) 

(59) 

(60) 

( 61 ) 


E 


E 


dfjz) df{w) 
dzi dwi 
dfiz) df{w) 


dzi 


dwj 


dEifiz)f{w)) 

dwi 

d^E{f{z)f{w)) 

dzidwi 

d^E{fiz)f{w)) 

dzidwj 


= ZiC 


= (1 + ZiWi)e^"^, and 


= ZjWie^"^ for i j . 


We now use (58-61) to compute that the only non-zero covariances in Cn are 


(62) 

(63) 

(64) 

(65) 

( 66 ) 

(67) 

( 68 ) 

(69) 


a 

6 


7 := 

p := 
ly := 

ri := 


E{f{x)f{x)) = E{f{y)f{y)) = e ^ , 

( dfiy)df{y) \ ^^( dfjy) df{y) \ g 

V dy. dy, ) V dy, dy, ) ^ ’ 

^ ( dfjx) df{x) \^ ^ ^ ^ dfjy) dfiy) \^ ^ g# 


dXn dXn J y ^yn dyn J 

E{f{x)f{y)) = e-"-, 


dfjx) 

dxn 


2 ® 


f df{x) df{y)\ ^ , 

E —---- = e 4 for z ^ n, and 

\ ^Ui ) 


uXn ^Vn J 


This proves that C„ also has the structure stated in Lemma 17. 


□ 


To apply the Kac-Rice formula to compute Kn^d{x, y) and Kn{x, y) for the values of x and y given in (42) we 
will need to compute det C„ d, det Cn, Jln.d, and the diagonalizations of Lin.d and Jl„. Here, Jln,d and are 
the matrices obtained from and C“^, respectively, by deleting all of the rows and columns whose indices are 
congruent 1 modulo n -I- 1, as in the Kac-Rice formula. 

By Lemma 17 we can do all of these calculations in terms of the generic form C given in (43) and then substitute 
in the values of a through r from (50-57) and (62-69) accordingly. 

The determinant of C is 
(70) 

det(C) = (/3^ — 77 ^) ^ ^ — ar — 6'^ — 2 S v + J fi — fiT — (a^ + ar — + 25v — ^ y —yr — 1 /^) . 

Recall that C = diagn(C) where C is described in (43) and (44). Applying a suitable permutation to the rows and 
columns of C, one obtains a block matrix with one 4x4 block and n — 1 copies of the same 2x2 block. Because 

of this, will have the same block structure and it can readily be computed to be = 


D+ 


E+ 


where Dh 
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and E± are the following (n + 1) x (n + 1) matrices; 


D+ = 


2 2 r2 of 2 

a "7 —CKT —d ” 7—2 o i/* r —7 iv" 
A 

0 


p 


a 6 v T — d^—S fj, r+d — 'y fi u n 

+ A 


E± = 


P 2_^2 


± 


Q: (5 T + Q: 7 1/ — S -y fl — U T—l'^ 


n a (5 7 + 0 ! T—(5^ —(5r+(5 — 7 /i 1 / 

U + ^ 


q:"^ 7 —a 5'^ —a iv''^+2 <5 iv -—7 /i"^ 
A 


T 


Q^T + Q:7I/' + ^^t/ —^7// —//t/T— 




n Q^r+2 a <5 — ^ 

^ A 


where 

(71) 

A = 0 ^ 7 ^ — c? — 2 o ^^7 — 4Qi^//r — 2 a 7 i/^H-^'^ + 2 r — 2 + 4 ^ 7 /i z/ — 7 ^/^^ + + 2 z/^r + z/^. 


If rj is obtained from C ^ by deleting all of the rows and columns whose indices are congruent 1 modulo n + 1, 
as in the Kac-Rice formula, we have 


Lemma 18. 17 = ^ 


(73) 


( O 1 with Cl = 

^^ 1,1 Hi ,2 

, where: 





ll 2 ,l ^ 2 , 2 . 




2 = diag 


/3 

2 

a 7 

) 

— a 6'^ — a h''^ + 2Sfiiy — 


to 

1 

772 ’ ■ ■ ■ ’ /32 _ 

9 ’ 

Tf^ 

A 

5 


V 

n —1 times 


) 



/ 





1 = diag 


V 

V 

a^T + 2a6v — S^fi — jjF’r — 

/32 

- 772 ’ ■ • ■ ’ 732 _ ^2 ’ 

A 



V 

n —1 times 



/ 


We notice that there exists a permutation matrix Q such that 

/ 


(74) 


M := Q'''J7Q = diag 


Ml,..., Ml, M2, ..., Ml,..., Ml, M2 


n— 1 times 


n— 1 times 


n times 



r p 

7? 


Q;^ 7 —q;l'^+2(5//L'— 7 /i^ 

cP" T-\-^OL 8 v — 8^ yu—p?" T— 

where Mi = 

p'^-v'^ 

V 


and M 2 = 

Q;^r +2 a. 6 — 

cP^ — oi 5^ — ol 1/^2 (5 7 /j,^ 


1^2 —rj2 

132 —T}^ 


L A 

A J 
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Lemma 19. We can orthogonally diagonalize with QP, where Q is the permutation matrix described in the 


previous paragraph and P = diag „2 


1 



i 

2 

2 

2 

V2 

2 


, obtaining 


(75) A := QPTr2QP == P^MP = diag 


Ai, A2,..., Ai, A2, As, At, ..., Ai, A2,..., Ai, A2, As, At 


n— 1 times 


n—1 times 


where 

(76) 


Ai — 


1 


P-v 

and At = 


A 2 — 


1 


As — 


a + fj. 


P + r]’ ay — ar — — 2i/S + yp, — p,T — 

a — pL 


ay + ar — 5'^ -\- 2v 5 — y p, — p,T — v'^ 


We now begin substituting in the values of a through r given in (50-57) and (62-69) into the results we have 
obtained for the C in order to derive the results we need for C„ ^ and C„. 

Lemma 20. For all t > 0, C„ is positive definite. For d > 3 and sufficiently small f > 0, Cn,d is positive definite. 

Proof. It is a general fact from probability theory that the covariance matrix of a random vector is positive semi- 
definite. For Cn, equation (70) becomes 


det(C„) = (e2* -e 2 ‘j -e 2 * /g 2 


±4.^ L 4.^ O 

2 * -e" 2 * 


y. 


which is positive for all f > 0 . 
For Cn,d we have 


aeTC..,) = ^ 


□ 


which is positive for d > 3 and f > 0 sufficiently small. 

Lemma 21. and fin are orthogonally diagonalized by QP where Q and P are the 2n^ x 2rP matrices described 

in Lemma 19. The eigenvalues \d,i^ Ad, 2 , Ad, 3 , Ad ,3 offln,d satisfy 

(77) 

= \/^^< + 0 ((=) A-f = ^/M + o«) x-f = yd(rrT)t x-f = ,pLl^±Ae + 0{ty 

The eigenvalues Ai, A 2 , A 3 , At of fin equal 

(78) 

.2 .2 .2 .2 

-1 


A, =(et 


A 2 = e 4 


/ li _ii\ 
(^e 4 -p e 4 j 


2 N -1 


As = 


6 4+6 4 


^2 _|_ e 2 — e 


and At = 


6 4—6 4 


PZ +z 

— + 6^ — 6 5- 


and satisfy 


(79) 


= ^ + 0{t^) A-'/" = ^ + 0(f2) A-'/" = f + O (f3) xy^/^ = + O (t^). 




Proof. This follows from Lemma 19 and by substituting the values of a through t from (50-57) and (62-69) into (76). 
(The asymptotics in (77) determined using the Maple computer algebra system [1]. However, they are simple enough 
that one can check them by hand.) □ 

We will need the following calculation in Section 5; 

Lemma 22. We have 
(80) 

(Ai,dA2,d) (A3,dAT,d) ^ 


xj det Cd,r; 


= d-5f-" + 0 (f-”+ 2 ) and = i-" + o(i-"+ 2 ). 

Vd6t Cn 
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Proof. Using (70) and Lemma 19 for the generic form of the covariance matrix C we have 

(A3A4)-^’ 


(81) 




Vdet C„ 

The result then follows by substituting (50-57) and (62-69) and doing an expansion. 


□ 


5. Proof of Part (1) from Theorem 9: short-range asymptotics 

We apply the Kac-Rice formula to the covariance matrices C„ ^ and C„ and the submatrices Vln^d and r2„ of their 
inverses, as computed in Section 4. It applies because, by Lemma 20, C is positive definite for all f > 0 and Cn,d is 
positive definite for all d > 3 and sufficiently small f > 0. The proof will be the nearly same for each, so we will work 
with K.n,d{t) and then explain what change needs to be made for ICn{t) at the very end of the section. 

We apply the diagonalization of A = A„ ^ = (QP)^r2„ ^(QP) from (74) and (75) to the Kac-Rice formula (41) 
to obtain 

1 


(82) 


K-n^d (t) = 


( 2 ^)n(n+l) p^(^x)pd{y) ^/det Cn,a 


j I det ^ (x)|| det J 7 (t)| e 


where r := [ ri T2 ... t„ where n = [ np ri ,2 ■ • ■ n^n ] for each 1 < i < n. In 

these new variables, ^ and rj become the new matrices ^ (t) and r] (x), whose entries are defined by 

^ V 2 

(83) (r) = — {-npj-i -I- Tj, 2 i) and (r) = — in^ 2 j-i + n^ 2 j) for i,j < n. 

The reason for diagonalizing r2„,d was to change the exponent into a form conducive to forming n sets of 2n- 

n 

dimensional spherical coordinates so that (Ar, r) becomes ^ rj.. 


Let !(;:=[ ri r 2 


k=l 

61 02 ... On ], where 0 i = 


0i,2 




(84) 


K.d'^i f n sin 6 »i_fe j {cos0- 


Ti,i = 


^2,d''’i 


\ ^ r ■ 

'^34''^ 


n—1 


\,d'^i 1 n sin6>i,fc ) (cos0i_„) 


fc=i 

n-l+i N 

n sin6>i,fc 

k—1 j 

'2n-l 

n sindi^fc 

fc=i 


cos 6, 




if j is odd and j ^ 2n — 1, 
if j = 2n, 

) if j is even and j ^ 2 n, and 
if j = 2 n — 1 . 


Thus 


IL 

, dr becomes (Ai^dA 2 ,d)~^"^”~^^ dfj,(^§,2n-i^«dr , where S' 


2n-l 


denotes the unit 


1=1 


sphere in and d^(g 2 n-i)n denotes the product measure on (S^" obtained from the standard spherical measure 

_ 1 

d/rg2n-i on S^”“^. Let cjiij be the trigonometric product in tij so that tij = After this variable change, 

we see that 


V2 




(85) e,j (r',0) = —r, 
where m = 1 when j ^ n and m = 3 when j = n. 


and (r, 0) = 


-1 + A, 


2 i-i 




Thus, in the new spherical coordinates, we have; 
( 86 ) 


^ {^l,d^2,d) (A 3 dA4 d) /" I J i /^^ II J i ! /^^ I 2 ?, ’’^TT 2ri-l4 4 

t t 1 lAi / Idet^(r,0) ||detr,(r,0) |e '==1 d/i(g 2 „-i)ndr. 

(27r) Vd(a;)/Od(y)Vdet i=i 

Using the asymptotic behavior of , ^ 2 ]/'^ ^ '^3 d^^’ expressed in (77), we notice that each the ele¬ 

ments of the nth column of each determinant vanishes linearly with t. Therefore, we factor out t from each column to 
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prevent these columns from vanishing in the limit as t goes to 0. Also note that each element of row i in both ^ and 
Tj are linear with ri. Thus, we let ^ (0, t) and fj {9, t) denote the resulting matrices when t is factored from the nth 
column and is factored from each row of each matrix, ^ and rj, respectively. Using Fubini’s Theorem, we can now 
split the integral from ( 86 ) into an integral over the radii and an integral over the angles; 


(87) 


( 88 ) 


ICn,d (t) = 


i^l,d^2,d) (A3,dA4^d) 


(2^)n(n-El) p^(^x)pd{y) ^/det Cn^d 


( 


n 


, 2 n+l 

I 


e *"=1 dr 


t J I det I ( 0 , t) II det f) ( 0 , f) I d^(§ 2 n-i)Ti 

y(S2n-l)n 

Using the definition of the Gamma function, (87) simplifies to: 

^ 2 "'r(n + l)” W (A 3 .dA 4 .d)”^” 


(89) 


K.n,d (f) = 


pd{x)pd{y)^ 

( 


y v/det Cfi.d 

det ^ (0,f) II det f) (0, f) I (i/i(s 2 n-i) 


r 


By (80) we have that 


(90) 


(Al,dA 2 ,d) ('^3,dM,d) 


= d-^ + o . 


\J det Cri,d 

Meanwhile, by (77), each entry of ^ and i) is of the form; constant (potentially 0) plus 0 {t). Therefore, 


(91) 

where 

(92) 


det ^ (0, t) 11 det 77 (0, t) \ d/r(s2„-i)n = Dn + O {t), 


(§2n-l)" 


Dn = lim 

t —^0 


det i ( 0 , t) 11 det f) ( 0 , t) \ dp(^s■2,^-l)n. 


(S2n~l)r. 


Pd{x) = pd{y) =-K 


d® + 0 (t 2 ). 


Finally, we also have from ( 6 ) that 

(93) 

Therefore, 

( 2"^7r"+^r (r?-I-1 

(94) where ^ 

We now compute the constant From Equations (85), (77) and (79): 

(95) linr|ij- (0,f) = lim (0,f) = \fd for j < n, and 


Dn 


(96) 


lirn ( 0 , t) = lim fpp ( 0 , t) = 


t —^0 


t —^0 


d(d- 1) 


6 i, 2 n-i forj = n. 


d{d-l) ; 


is factored 


Let (1 (0) be the resulting matrix when \/d is factored out of the first through n — 1-st columns and 
out of the nth column of lim ^ (0, t). If we do the same process with lim p, (0, t), we obtain the same result with the 

t —>0 t —>0 

sign changed in the nth column. Therefore, 

d"(d-l) f 
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From (84), we notice that each entry of row i of (0) contains a factor of sin 9ij. Thus, we can take this factor 

j=i 

out of each row of the matrix, removing any dependence of the determinant on 9ij for j < n. Let v {6) denote the 
matrix that remains after removing these factors. 

We can then split the integral into two, an integral over B", where B := [0, tt], corresponding to Oij for j < n, and 
an integral over corresponding to O^ j for j > n: 

(98) D„ = /n dLebe" J | det iv ( 0 ) . 


Here we have used that d/r§ 2 n-i = sin 0^" ^ ^dLehs . The former integral in this product can be calcu- 

i=i 

lated recursively with integration by parts to be; 


(99) 


(sin 0 ij)^”’'‘^ ^ dLebB" = ( H ( / (sin 6 ») 




i,j^n 


Vi=i 



= 7rrtT2LtJ. 


2n!! 


The calculation of the latter integral follows from Proposition 30 of Appendix B; 


( 100 ) 




1 ). 


r(^) 


r(^) lr(H±4) 


r (n + 1) 7r"2 

rm" 


Therefore, 

( 101 ) 

' 2"V”+ir(n+ 1)" 


^n,tt 


( 102 ) 

Thus, 


n(n+l) p ^ n+1 ^ ^ 


,(27r) 


D„ = 


2^\^+^r{n + iy 


.(27r) 


n(n+l) p /'n+ 1^2 ^1 


d^id-l) \ y r(n + l)^ 

2nUj 


d-i\ v^r(i^) 

2 r(^) 


= a 


n,d’ 


d 1 ^ V^r( + ^ 


ICnM (i) - ( J 2 p 

as stated in Part (1) from Theorem 9. 

The only differences when computing K.n{t) instead of ICn^d{t) 

( 1 ) = \/2 + 0{t^) instead of v/M + O(f^) and =t + 0{t^) rather than + 0{t^), 

( 2 ) the factor of in (80) is missing, and 

(3) the factors of c? 2 ” are missing from the expression for the density of the zeros. 

One can readily check that this results in the factor of being removed from the constant: 

f-y (n+2\ 


□ (Part (1) of Theorem 9). 


6 . Proof of Part (2) from Theorem 9: long-range asymptotics 

It will be convenient to apply the Kac-Rice formulae to the ensemble g given in (14), which has the same zeros 
as the Isom(]R")-invariant ensemble f. Let Cn,g denote the covariance matrix applied to random vector (40) for this 
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ensemble. Recall that x and y are given by (42). The following covariances can be computed from those in (62-69) 
and the product rule; 


(103) 

(104) 

(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 


E{g{x)g{x)) 
dg{x) 


E g{x) 


E 


E 


dxj 
dgjx) dg{x) 
dxj dxj 
dgjx) dg{x) 
dxi dxj 


/ dgjx) dgiy) \ 
dx, dy, ) 


E 


= E {g{y)g{y)) = I, 

e( ^ 9{y) \ ^ ^ 

V dyj dyj ) 

^ (dg{y) dg[y)\ ^ . 


E{g{x)g{y)) = e = 


-o\\x-y\\ _ .- 


= e 2 




— c^-oW^-yw 


— Q-2\\^-y\\ 


(^1 - {Xi - yif'^ 


i.Xj -yj) = 


0 


if j ^ 


—te 2 if j = n. 
6 2 if i ^ n 


dgjx) dg{y) \ 
dxj dyk ) 


-e 2 


e 2 (1 — if i = n 


and 


Remark that Cn,g has the structure asserted in (17) and that 


, ^ / , 2 \ "("-!) / 
det(C„,g) = - e * j 


1 -be- 


e - e-‘ 




1 — e 


- 


> 0 


for all f > 0 , so that C„ g is positive definite. 
The proof will rely upon two facts: 


(111) (detC„,g) ^ = 1-f 0(f^e ‘ ) and 

(112) ||I-n„,g||oO =0(fV‘'), 


where 11 | |oo denotes the maximum entry of the matrix. The former can be obtained from expression (70). The latter 
follows from the calculations above and Lemma 18 expressing r2„ g in terms of the entries of C„ g. 

The covariance matrix for random vector (38) is the identity, by (103-106) above. Thus, the Kac-Rice Formula for 
the density of zeroes of gix) gives 


1 = 


p{x)p{y) ^ 1 

p{x)p{y) p{x)p{y) 
1 


1 


(271-) 


»(rt + l) 




1 


(27r) 




|detJ 7 |e 2 (n.u)(^^ 


(27r)"'(”+i) p{x)p{y) 


J I detail detryje 2 


(113) 
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Since C^n.g is positive definite, the Kac-Rice formula for two-point correlations (41) and Equation (113) give 

J I detail - 1 


\}Cn{t) - 1 | = 

(114) 


^ (27r)”^" p{x)p{y) y^det 
1 


R2" 


(115) 


(116) 


< 


( 27 r)"-("+i)p(a;)p(y)y^det 

J I detail _ y'det C„,g J | det||| detr/|e-5(“'“)d. 

R2n2 

dct Cn,g 

(27r)"'("+i) p(a;)p(y) 


u 


R2" 


J I detail detr/je _ J | det||| detryje 


R2" 


1 - (del Cn,g) 2 


(27r)"("+l)p(a;)p(y) 


J |det|||detr 7 |e 


R2" 


Equation (116) is by (111). 

From Lemma 13 Part 2, using A = I and B = ^n,g, we have 


(117) 


J I detail detr7|e-5(“".s“-“) - J | det||| detr/|e-5(i“’“)dit 


o(\\i-n^,g\\U^)=o{te-^). 


Therefore, |/C(f) — 1| = O (fe *2 


so we obtain the desired result. 


□ (Part (2) of Theorem 9). 


7. Proof of Theorem 3: Universality 

This section is devoted to a proof of Theorem 3. It will be divided into three parts: 1) Reduction to a local version, 
2) Statement and proofs of two lemmas, and 3) Proof of the local version. 


7.1. Reduction to a local version of Theorem 3. Let the homogeneous coordinates on RP^ be denoted [Zi,..., Zk+i]. 
After applying a suitable isometry from SO{k + 1) we can assume that p = [0 : • • • : 0 : 1], allowing us to work in 
the affine (local) coordinates 


(118) 


^1 

Zk+i ’ 


Zk 


Zk 

Zk+i 


having [0 : • • • : 0 : 1] G RP^ as their origin. 

In these local coordinates, the tangent space Tp{M) becomes an n-dimensional linear subspace of R^ containing 
the two points x ^ y. We can now rotate about 0 in these local coordinates by an element of SO{k) (corresponding 
to an element of SO{k + 1) that fixes [0 : • • • : 0 : 1]) allowing us to assume that 

(1) TpM = span(ei,..., e„), where ei,..., are the standard basis vectors on R^', and 

(2) a: = (0,..., 0, s, t) and y — (0,..., 0, u) in the local coordinates (zi,..., Zn) on TpM = R", 

where {s,t) ^ (0,u). Since the ensemble (3) on RP^' is invariant under elements of SO{k + 1) and since we 
have rotated the submanifold M and the points x and y under the same composition of elements of S'0(fc + 1), the 
correlation function remains the same. 

By our choice 1, above, M is locally expressed as a graph of a function xjj : R" R^“" that satisfies 


(119) ■0(0) = 0 and 790(0) = 0. 

The orthogonal projection proj^ : TpM = R" —>■ M is given by projp(m) = ((m), 0(nj)) for any w G R". 
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In affine coordinates (118), the SOik + l)-invariant polynomials are 

( 120 ) 


fdiz) = 

\a\<d 



acyZ 


where 


dl 


{d- |a|)! n cti' 

i=l 


and the are iid on the standard normal distribution A/^ (0,1). 
The correlation between zeros 


K.. 


n,d,M I PrOJp 


X 

yfd 




is the same as the correlation between zeros for the pull-back of this ensemble to the tangent space TpM = 
under proj^, which is given by systems of n functions chosen iid of the form 


C 


( 121 ) 




Vd 


■= fd 


Vd 




Vd) 


This follows because one need not use round balls in the definition (8) of the correlation function-any sequence of 
neighborhoods that is sufficiently nice for computing a Radon-Nikodym derivative suffices (see Remark 12). If one 


uses round balls Ns (projp ) and iV^ (^projp (;^) ) i 


in the definition of K, 


NM (projp(;^):Projp(^)), 


then their preimages under the mapping projp will be suitable neighborhoods for defining the correlation function 
for the pull-back (121). 

Thus, we have reduced the statement of Theorem 3 to: 

Theorem 23 (Local version of Theorem 3). Let = Kd.tp denote the correlation function for systems of n 

functions chosen iid of the fonn (121) and let denote the correlation function for the Isom(K")-invananf sys¬ 
tem (4). 

For any s, f, u S K., i/a; = (0,..., 0, s, t) G K.", y = (0,..., 0, it) G R", and x y, then 

with the constant implicit in the O notation depending uniformly on compact subsets ofMf x R \ {(s, f) = (0, it)}. 

We will need the following more detailed notation in the next two subsections. Letip{x) = {'fi{x),... ,^fe_„(a;)). 

(.1 and 7 G Z’fT 

fd 


K, 


n,d.'ip 


If we write a = (/3,7) with /3 G and 7 G ”, then (121) becomes 


( 122 ) 




\Vd) 


E ^ 


iP,-r) 




|/3|-|-|7|<ct 


\VdJ 


where b 


{13,f) 





and 


if 

d 

(3 7, 


s/d J 


dl 


n k—n 

{d- |/3| - I 7 I)! n A! n 7*! 

i=l i=l 


As before, the coefficients are iid on the standard normal distribution Af (0,1). 

7.2. Two lemmas. 

Lemma 24. For any x = {0,0,..., s,t) and y = (0, 0,..., 0, u) in R” with x y we have: 

(1) The covariance matrix C corresponding to random vector (38) from the Kac-Rice formula for density (39) 
applied to the Isom(M")-tnvanflnf ensemble (4) (at x or y) is positive definite. The submatrix fl of C~^ 
defined in (39) is also positive definite. 

(2) The covariance matrix C corresponding to random vector (40) from the Two Point Kac-Rice formula (41) 
applied to the Isom(R”)-tnvananf ensemble (4) is positive definite. The submatrix fl ofC~^ defined in (41) 
is also positive definite. 

Proof. We give the proof of Part 2 leaving the necessary modifications for Part 1 to the reader. 

It is a general fact from probability theory that the covariance matrix of a random vector is positive semi-definite. 

Thus, it will be sufficient for us to check that det(C) > 0. We substitute a: = (0,0,..., s, f) and y = (0,0,..., 0, it) 
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into the covariances computed in Equations (58-61) from the proof of Lemma 17, obtaining that C is of the form 
A B^] 

and A, B, and D are the following (n -f 1) x (n -I- 1) matrices: 


diag„(C), where C = 


B D 



1 

1—‘ 

o 

1 


'10.. 

1- 

o 


0 1 ■ 

0 


0 1 

0 

(123) A = 

s 

_t 0 . 

1 -1- st 

st 1 + 

, B = e*“ 

s 

f 0 .. 

1 su 

0 1 + tu 


and 


(124) 


D = 


1 0 
0 1 


0 

u 0 


u 

0 


0 1 +u^ 


After applying a suitable permutation to the rows and columns, C becomes a block-diagonal matrix with n — 2 copies 
of 


gtu 


and one copy of 


(125) 



se + 

fe®'+*' 

e‘“ 

se‘“ 

fe‘“ 


(1 -1- S^) 6®"'+*'’ 


0 

e*“ 

0 

2 1 .2 
fe® +* 

ste^ +‘ 

(l+t2) 

ue*" 

sue*'^ 

(1 -1- tu) e' 

e*“ 

0 

ue‘“ 

e“^ 

0 


se*“ 

e*“ 

sue*"^ 

0 


0 

fe*“ 

0 

(1-f fu)e*“ 


0 

(l -1- e' 


The former has determinant e^‘“(e® — 1), which is positive since r = s'^ + (t — > 0, by our hypothesis 

that X ^ y. The latter has determinant equal to 

Without the exponential prefactor, this equals (e’’ — ’’) (r^ -f 3) -f — 1, which one can also check is positive for 

all r > 0. 

Since C is positive definite, so is C“^. After applying a suitable permutation to the rows and columns of C~^ the 
X n? principal minor is 17, which is therefore positive definite. □ 


Lemma 25. For any x = {0,0,..., s,t) and y = (0, 0,..., 0, m) in K." with x ^ y we have: 

(1) Let Cd.tp and C be the covariance matrix for vector (38) applied to the systems hd,ip given in (121) and f 
given in (4), respectively, (at either x or y) and let Fid,ip tind FI denote the submatrices o/C^^ and 
defined in the Kac-Rice formula for density (39). Then, 


C 




c + o 



and 


Fld,tp — 17-1-0 


where the constants implicit in the notation depends uniformly on compact subsets (if we are working at 

a:) or K. (if we are working at y). 
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(2) Let Cd,xi; o.nd C be the covariance matrix for vector (40) applied to the systems given in (121) and 
f given in (4), respectively, and let fld,tf) and H denote the submatrices o/C^^ and defined in the 
Kac-Rice formula for density (41). Then, 


Cd,^ = C + Oi- 


and 


^d,tp — + O I — 


with the constant implicit in the O notation depending uniformly on compact subsets o/K^ X K \ {(s, t) = 
{0,u)}. 

Proof We will prove Part 2, leaving the necessary (simple) modifications for Part 1 to the reader. 

We will only use the assumption that x = (0,0,. .., s, f) and y = (0,0,..., 0, u) with a; ^ y in the last three 
lines of the proof, to obtain the estimate relating fid.-i/i to 17. Until then, x and y will denote any two points of R". In 
particular, the estimate relating Cd,ii, to C holds for any x,y £ R". 

We will first show that 

(126) Cd,rp = Cd,o + 0 



Because tp is C^, there exists a constant ^ > 0 independent of d such that for any multi-indices (3 € Z" and 7 € Z^ " 
we have; 


(127) 

(128) 

(129) 


a 

dxi 




Both the constant A and the multiplicative constants (implicit in the O notation) depend uniformly on x within compact 
subsets of R”. 

We can now prove that for any x,y G R" 


(130) 



ItI^i 






where we can have x = y and/or i = j. The constant implicit in the big-O notation depends uniformly on compact 
subsets of R” x R". 
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The proofs of each of these are essentially the same, so we’ll prove (132), leaving the proofs of (130) and (131) for 
the reader. We have 


(133) 


d ( / X 


dxi 




Ip 


VdJ V K'/d 


d 

9yj 




. ((JL)'’) (JjDYjlI (jlY) U (X) 


+ 


dxi yxVdJ 

p 


p'^ 


Ip 


VdJJ dyj y\y/dj 

f y Y ^ 

Vd)) \Vd) dU] 


VdJ 


Ip 


y 

y/d 


+ 


+ 


Vd 


VdJ dxi 

X Y ^ 
Vd) 


%p 


Vd) J ) dyj \ \VdJ 

(jlY — 

\Vd) dyj 


Vd 


Ip 


< B 


y 

Vd 

^2\ I/sI+sI-yI 

tJ 


Here, A, i? > 0 are given by (127-129) and are independent of |/3|, I7I and d, but depend uniformly on x and y within 
compact subsets of R” x R". In particular. 


E 


l^l + |-rl<<^, 
|7|>1 


d 




P 


Ip 


PjJ dxi \ \VdJ \ \Vd)) I dyj \\VdJ \ \Vd) 


tX g / / „. \ P / / „. w-y' 


fxpf 


^ E 

l3l + lTl<d, 

IyI^i 

BA^ 


[dlj^\d) 


^2^l/3|+2|'rl ^^2 

x: 


\m+\~i\<d, 


d \ 


/3 7/ \ d J 


< 


E 


\P\ + \-t\<d 


filJ \ d) 


^2nI/3| + N ^^2/ a 2 


< 


1 + k- 


< 


C 


for some 


C > 0, since limd_>oo (l + 


= e 


kA^ 


The estimate (126) follows immediately. For example, 

dhd^^jx) dhd,-^{y) \ _ 

dxi dyj J 




\P\<d 


IyI>i 


= E 


dhd^ojx) dhd^o(y) \ ^ q Pl_ 
dxi dyj J \d 


We will now show that 
(134) 


Cd,o = C + 0[-] . 


It follows from a calculation analogous to that from Lemma 17, but with a rescaling by IjVd and the fact that 




= e^ + 0 


with the constant depending uniformly on a; € R. Rather than including the computation for each of the eight different 
types of expectations, we simply list one of the more complicated ones here: 
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E 


/ dhd.ojx) dhd^o{y) \ 
V dxi dyi ) 


E 


dfjx) df{y) \ 
dx^ dyi ) 


x y Xiyi{d—1)\ / x-y\d--^ 

(1 + Xiyi)^^"^^. 


Combining (126) and (134) we conclude that = C + O (i). 


If a; = (0,..., 0, s, t) and y = (0,..., 0, u), with x ^ y. Part 2 of Lemma 24 gives that C is positive definite. 
Then, there is a neighborhood of C in the space of all 2n^ x 2n^ matrices on which taking the inverse is a differentiable 
map. Therefore, + O (^) and hence r2d,i/> = ^ 


0{^d) 


□ 


7.3. Proof of the local universality theorem. Proof of Theorem 23. Throughout the proof we will use that we have 
normalized so that x = (0,..., 0, s, f) and y = (0,..., 0, u) with x ^ y. 

Part 2 of Lemma 24 gives that the covariance matrix C for random vector (40) applied to ensemble (4) is positive 
definite. Therefore, Part 2 of Lemma 25 gives that covariance matrix for random vector (40) applied to ensem¬ 
ble (121) will also be positive definite for d sufficiently large. We can therefore apply the Kac-Rice formula (41) to 
ensemble (121) obtaining 


(135) 


Kd,- 4 , {x,y) = 


(2^)n(n+i) pJ^,^(^x)pd,■,|,{y)^/detC^ 




R2" 


We will show that (135) differs by O from the result obtained when applying the Kac-Rice Formula (41) to the 
Isom(IR")-invariant ensemble (4). 

Let us first consider the prefactor from (135). To show that pd,tpix) = p{x) + O we apply Lemma 13 to 

the Kac-Rice formula for density (39). This follows because the matrix ft in (39) is positive definite, by Part 1 of 
Lemma 24, and because fld.ii, = ft + O (^), by Part 1 of Lemma 25. The same estimate holds at y. Meanwhile, 
since Cd,ii, = C + O (^), with C positive definite, it follows immediately that 

1 

V det Cd 

We now consider the integral in (135). Part 2 of Lemma 24 gives that the matrix ft used when applying the Kac- 
Rice formula (41) to the Isom(R")-invariant ensemble (4) is positive definite. Meanwhile, Part 2 of Lemma 25 gives 
that ftd^ij, = r2 -f O (i). Thus, Lemma 13 gives that the integral from (135) differs by O from the integral in 
(41), when (41) is applied to (4). Each of the lemmas used asserts that the constants depend uniformly on compact 
subsets of X K. \ {(s, f) = (0, u)}. □ 


VditC 


+ 0 - 


Since we reduced the statement of Theorem 3 to Theorem 23 in Subsection 7.1, we have also finished the proof of 
Theorem 3. 


8. Finite degree restrictions to submanieolds M c MP^ depend on the geometry 

We present a simple example illustrating that the constant from the leading term in the correlation function can 
depend on the geometry of a submanifolds M C RP^ if the degree d of the ensemble is finite. Thus, it is not possible 
to prove a universal formula for the short-range asymptotics at finite degree. 

We consider the ensemble S'0(3)-invariant polynomials of degree 3 because for degrees 1 and 2 the covariance 
matrix for random vector (40) is not positive definite. Restricted to a system of affine coordinates (x, y) ^ [x : y : 1] 
each such polynomial has the form 

F^^x, y) = 00,0 + y/iaifiX + v^oo.iy -I- V^aipxy -I- V^a 2 ,ox^ + VSao^ 2 y^ 

+VSa2,ix'^y + v^ai,2Xi/^ -I- 03 ,qx^ -|- ao,3y^, 

where each of the coefficients is chosen iid with respect to the standard normal distribution, Af (0,1). 
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Let M C be given by ?/ = bx^. As in the previous section, we parameterize M by the x coordinate, in this case 
forming the one-variable ensemble of random polynomials that depend on 6 as a parameter; 

(136) Hs,aix) = 00,0 + VSaifiX + Viao,ibx^ + VQai^ibx^ + V^a2,ox^ + Viao,2b^x‘^ 

+y/3a2,ibx^ + •\/3ai,26^x® + Qs^qx^ + ao^sb^x^, 

In order for our results from Sections 4 and 5 to apply here, we multiply by the prefactor: := 

We will use the Kac-Rice formula (41) to show that the value of b affects the constant term in the short-range 
asymptotics for the correlation function K{x, y) with x = —^ and y = |. 

An easy calculation using the Kac-Rice formula (39) for density of zeros gives that 


p{x) = p{y) 


Vs 


oin 


where the constant term agrees with the result ( 6 ) that is stated in the introduction for the S'O(2)-invariant polynomials 
of degree 3. 

The covariance matrix for random vector (40) applied to the ensemble with x = — | and y = | is 


(137) 


C = 


a d p —V 

S "f ly T 

p u a —6 

—V T —6 7 


where 


-f 4f2 + 16)' 


<5 = - 


3f {kH"^ + 2) {kH^ + 4f2 + 16)' 


7 = 


4096 ’ 1024 

{SW + 12f2 + 48) {SkH^ + 13fc2f4 + 16^2^2 X2f2 + 16) 


k- = 


(fc2f4 4^2 X0)'^ 


256 

and 


4096 

3t {k^t^ — 2) {k^V — 4f2 4- 16)^ 

’ “ 1024 ^ 

{ikH^ - 12f + 48) {SkH^ - UkV* + 16kH^ + Uf - 16) 

256 ' 


This matrix has exactly the same structure as that from Sections 4 and 5, in the case that the dimension is 1. Therefore, 
the submatrix of C~^ from the Kac-Rice formula (41) is diagonalized in precisely the same fashion, with the 
eigenvalues satisfying 


Ag = V 6(1 + ^ A4 + 362 ^2 _|_ o{V). 

(We call them A 3 and A 4 in order to be consistent with the previous sections.) We also have 

detC= (54 6^-^63 62-4 9) V + 0{t^°) . 

The calculation of the short-range asymptotics done in Section 5 applies here, with the minor modifications to 
Ag , A 4 , and det C listed above. One obtains 

Proposition 26. The correlation between zeros for ensemble (136) satisfies 

In particular, the leading term depends on the cur\’ature of M at 0. 

Question 27. In the general setting of M G RP^ how does the constant in the leading order asymptotics near p G M 
depend on the local geometry of M at p? 
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9. Universality in the complex setting 

We begin by adapting the Kac-Rice formulae (39) and (41) to the complex setting. As the modifications are nearly 
identical, we will discuss the formula for correlations, leaving the formula for density to the reader. 

Suppose that h = {hi, / 12 , ■ • ■ j : C" —?► C" is a Gaussian random analytic function with complex Gaussian 
coefficients. Let £ and rj be the n x n complex matrices whose rows are ^ 1 ,... and r/i,... r/n, respectively. Let 
"t* = [ Vi ^2 172 $,n ] ^, the vector formed by alternating the vectors and ? 7 i. 

Proposition 28. Suppose that the covariance matrix C = E{viVj) of the random vector (40) is positive definite. Then, 
the two-point correlation function for the zeroes of the system h is: 

(138) y) = "(I- 

where ft is the matrix of the elements ofC~^ left after removing the rows and columns with indices congruent to 1 
modulo n 4- 1, * denotes conjugate transpose, and (, ) denotes the Hermitian inner product. 

Our General Kac-Rice Formula (Proposition 10) can be easily adapted to the complex setting, from which Proposi¬ 
tion 28 follows easily. Alternatively, Proposition 28 follows from [8, Theorem 2.1] by using the suitable Gaussian 
density Dk{0, z) in their formula 32 and normalizing by the density at the two points x and y. 

With these modifications to the Kac-Rice formulae, the proof of Theorem 4 adapts nearly verbatim from the proof 
of Theorem 3. We list below the simple modifications that need to be checked; 

(1) The proof of Lemma 13 adapts easily to the integral expression in (138) and to the analogous formula for the 
densities p{x) and p{y). 

(2) The proof of Lemma 24 is easily adapted. More specifically, the determinant of the 6x6 block analogous to 
(125) equals; 

g6Re(tII) + (|(-.|)^ _ g2(|s|) = +2(|t-.|f ^ ^ ^ ^ ^3 (L|f+3 (|i-.|)^ _ 

which is positive for (s, t) f (0, it). 

(3) The proof of Lemma 25 applies after verifying that, when expressed in local coordinates, the covariance matrix 
for the rescaled SU{n l)-invariant polynomials differs from that of the Isom(C")-invariant ensemble by 
O (i). (These covariances are listed in [8, Sections 2.4 and 4].) 


Appendix A. Proof of Lemma 13 
We will need the following lemma to prove Lemma 13. 

Lemma 29. For any positive definite mv? x m.'nf matrix A, there exists a constant D > 0 such that, for R sufficiently 
large. 


(139) 


/ 


2=1 


£ e 


-DR 


(Here u is as in (23). 

Proof. First, we note that (Ait, it) > Aminll'*^! P- The left side of inequality (139) is bounded by 


(140) T 

using that U™ i I 


n 

det^l < n!'"| 


det|*|e 

itll™" . Let a ;= 2mn^ 


OO 


- 1 and b := 


(141) 






_ —6r^+ar 


dr = 


V^e- 


-erfc 


R 


R 
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The result then follows because erfc(x) := —j= / e ^ dt < e ^ 


□ 


Proof of Lemma 13. We first split the left side of the inequality into integrals with ||tt|| < R and ||m|| > R for 
some R: 


(142) 


(143) 


(144) 


]^|det|*|e = 






2=1 


2=1 


< 


/ m 

II 


|2x||<i? 




|2x| \>R 






For B sufficiently close to A, min (Bm, m)y/ 2 > (Art, u). Thus, Lemma 29 gives that the two latter summands of 
(142) are both bound by e~^^. 

We use Holder’s Inequality to bound the first summand in (142) with 
(145) 




|u||<_R 




L“(||«||<fl) 


= D5 




L~(||ix||<_R) 


Let L = A — B. Since ||e 2 “ — 1|| <a;asa; approaches 1, 


(146) 


,5(Lu,ll) _ 2^ 




< ll(L«:“)llL~(||«||<fl) < l|L|loo ll“ll||«||<fl < mn^R^\\L\\oo ■ 


Therefore, we have 


(147) 


J I detail detryle-^Ai^.^) - l) du 


< mn^D^R^ IIA — B| 


2e 


-DR 


The result follows if we set i? = IIA — Blloo^^^. 


□ 


Appendix B. Moments of the Volume of a Random Unit Parallelotope 

The derivation of the formula for the density (7) and the short-range asymptotics from Theorems 1 and 2 require 
the following formula: 

Proposition 30. Consider n random unit vectors in M" chosen independently with respect to spherical measure. The 
kth moment of the volume V of the parallelotope formed by these vectors is 




u_i r(i) 


_ ■ n£(^ 

r(^)y r(i) 


This result is proved in greater generality in [27] (see Equation 23); however, we will give a simple derivation of 
the formula. 


Proof The method for finding this moment involves fixing each of the vectors Vi that determine the parallelotope one 
at a time with respect to the parallelotope in described by the previous i—1 vectors, weighting each newly added 
vector based on the probability of obtaining a given height off of the previous (i — 1) —parallelotope. 
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Thus, after fixing the first vector, we see that the height of the second vector off of the first vector in terms of the 
angle 9 off of the normal to the first is cos {9) and the probability density of obtaining this height is 

^n-2g 


sm 


/siir 

0 


n-2 


6 de 


Yet, from the next vector on, the vector can vary in two directions, along two different spheres, one of dimension i—1 
and the other of dimension n — i — 1, in order to maintain the same height. Thus, for the zth vector, we have that the 
probability density of obtaining each height cos 9 for an angle 9 off the normal vector of the base given by the first 
i—1 vectors is 

sin"“*“^ 0cos*“^ 9 


/sin" 

0 


9 cos*“i 9 d9 


To express the fcth power of the volume, we multiply the fcth power of each of the heights together. Thus, 

A 


n—1 


(148) 


(149) 


E[v'^]= n 


/7r/2 

/ sin"-'-i6lcos*-i+'=6»d6» 

0 


V 


tt/2 

f Sin" 

0 


9d9 



'JAIL' 

rm. 


n— 1 


n— 1 


11 I r(ii^)r(i) 


k -p / n— 1 + 2 '' 


r(P 


□ 
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